RESEARCH  IN  SEISMOLOGY 
W .  Stauder,  et  a  1 
St.  Louis  University 


AD-784  958 


Prepared  for: 

Air  Force  Office  of  Scientific  Research 
Advanced  Research  Projects  Agency 

1974 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


AD784958 


I 


AfGSK  -  [fc  -  7  4 


i  j  r. 


Semi-Annual  Technical  Report 
October  1,  1973  -  April  30,  1974 


ARPA  Order  No.:  1827-5 
Program  Code:  3F10 
Contractor:  Saint  Louis  University 
Effective  Date  of  Contract:  1  March  1973 
Contract  Expiration  Date:  April  30,  1975 
Amount  of  Contract:  $103,340 
Contract  Number:  F44620-73-C-0042 
Principal  Investigator:  William  Stauder 
Phone:  314-535-3300-Ext .  540 

Program  Manager:  Same 
Short  Title:  Research  in  Seismology 


D  D  C 


Sponsored  by 


Advanced  Research  Projects  Agency 


ARPA  Order  No.  1827-5 


AIR  FORCE  OFFICE  Of  CC'EMTIFIC  RESEARCH  (AFSC) 
NOTICE  or  T  TAi  TO  roc 
1!  ;  I:.--.-:  end  is 

a.  j  .  .  V  I'  .  ALA  190-12  (7b). 

Dic'riLi.  . 

D.  V.  TAYLOR 
TccYriica;  l.i'c.  .  i"  .:r 


Reproduced  hv 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U  ^  Department  of  r <  mme-re 
Springfield  VA  27151 

ilk 


in 


REPORT  DOCUMENTATION  PACE  ; 

1  READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

'  <1  A  ,  ,  _  2.  GOVT  ACCESSION  NO 

-  7R -  <’4^1337 

13.  RECIPIENT'S  CATALOG  NUMBER 

4.  Tl  f  LE  (find  Subtitle) 

RESEARCH  IN  SEISMOLOGY 

5.  TYPE  OF  REPORT  6  PERIOD  COVERED 

Semi-annual  Technical 

L  October  73-30  April  74 

*•  PERFORMING  ORG.  REPORT  NUMBER 

W.  Stauder,  R.  Herrmann,  R.  Street, 

R.  Somayajulu 

••  CONTRACT  OR  GhF'IT  NUMBERfaJ 

F44620-73-C-0042 

».  PERFORMING  ORGANIZATION  NAME  AMO  ADDRESS 

Saint  Louis  University 

St.  Louis,  MO  63108 

10.  PROGRAM  ELEMENT.  PROJECT  TAS< 

ARFa  A  WORK  UNIT  NUMBERS 

ARP A  Order  1027-5 

Program  Code  3F10 

Program  Element  62701 E 

U.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Advanced  Research  Project  Agency/NMR 

1400  Wilson  Blvd. 

Arlington,  VA  22209 

12.  REPORT  DATE 

1974 

IS  NUMBER  OF  RACES  . 

4-tr 

MONITORING  AGENCY  NAME  A  ACORESS <•!/  F.H.IMI  Iron,  Contrail, nt  Oil, ro) 

Air  Force  Office  of  Scientific  Research/NP 
1400  Wilson  Blvd.  / 

Arlington,  VA  22209 

15.  SECURITY  CLASS,  (ol  Ihle  repot!) 

J  Unclassified 

IS«.  DECLASSIFICATION.' DOWNGRADING 
SCHEDULE 

i*.  Ulb  \  RI«U  MON  STATEMENT  (ot  thie  Report)  - — - 

Approved  for  public  release;  distribution  unlimited. 

IT.  OlSI  HI BU T ION  ST  ATEMENT  (el  the  ebetr eel  entered  In  Block  20.  II  dlllerenl  Iron,  Report) 

U.  supplementary  notes  "  '  - - - — - 

'*•  »uNo»  continue  on  reveree  eide  II  neceeeery  end  Identity  by  block  number) 

Source  parameters 

Focal  depth 

Seismic  Scaling  Law 

Central  U.S.  earthquakes 

Q  structure  at  continental  boundary 

Matching  theoretical  and  observed  Love  and  Rayleigh  soectra  allows 
determination  of  source  parameters  including  focal  depth.  Micro 
to  moderate  earthquakes  in  Mid-continent  scale  according  to  Aki's 
model  A,  and  evidence  a  two-fold  regional  trend,  north-south  and 
northwest.  The  0  structure  under  the  continental  boundary  of 

South  America  indicates  complex  structure  and  average  low  Q. 

SECURITY 


»ITV  £L  AS S!  FI r  A  T" ON  OFT h I S 


(H?ie/>  Dnrm  Entered) 


V  ..  i.u  (,i 


t  # 

1/ 


Research  In  Seismology 
Semi-Annual  Technical  Report 
Contract  F44620-73-C-0042 


Table  of  Contents 


Technical  Report  Summary 


1 


1.  A  Technique  for  Determining  Earthquake 
Source  Parameters  by  Using  Surface 
Wave  Amplitude  Spectra  2 


2.  The  Seismic  Scaling  Law  and  Source 

Mechanism  Characteristics  for  Small 

Magnitude  Events  in  the  Central 

United  States  11 


3.  The  Q  Structure  above  the  Subdue tion 

Zone  in  South  America  29 


Technical  Report  Summary 


This  research  concerns  magnitude  determination  and  the 
determination  of  source  parameters  in  a  region  character¬ 
istic  of  a  continental  interior.  While  the  ultimate  goal 
is  a  comparative  study  of  sources  of  like  magnitude  in 
Eurasian  source  regions,  the  proximate  purpose  has  been  a 
better  understanding,  first,  of  seismic  source  parameters 
in  the  mid-continent  United  States .  This  report  completes 
this  phase  of  the  investigation,  and  presents  results  in 
three  areas  of  inquiry. 

1.  A  technique  has  been  developed  whereby  theoretical 
Love  and  Rayleigh  amplitude  spectra  can  be  generated  for  an 
assumed  earth  model  for  a  given  choice  of  five  source  param¬ 
eters:  the  strike,  dip,  and  slip  angles,  focal  depth  and 
seismic  moment.  The  parameters  are  systematically  varied 

to  find  the  best  agreement  between  observed  and  theoretical 
spectra  for  the  entire  range  of  the  observed  spectra .  The 
method  has  been  applied  to  determine  the  source  parameters 
for  10  central  U.  S.  earthquakes  with  MQ  ='  io22  dyne-cm. 

2.  Other  studies  of  earthquakes  of  the  central  U.  S. 
have  resulted  in  a  seismic  scaling  law  for  earthquakes  in 
the  moment  range  10l8  through  1024  dyne-cm  which  resembles 
Aki's  revised  model  A  theoretical  law  but  displaced  toward 
shorter  periods.  Dependence  is  shown  on  the  period  at  which 


la 


mb  Is  determined  in  the  m^  vs  MQ  relation.  Focal  mechanisms 
of  38  of  these  events  show  a  complicated  regional  stress 
pattern  with  two  principal  trends:  a  north-south  trend  along 
which  high  angle  normal  faulting  changes  to  high  angle  re¬ 
verse  faulting,  and  a  northwest  trend  of  normal  faults. 

3.  The  spectral  ratio  of  pP/p  from  deep  focus  earth¬ 
quakes  in  South  America  yields  generally  low  values  of  § 

(Q  =  400)  for  the  region  above  the  focus,  with  systematic 
differences  for  segments  of  path  along  rays  to  western  U.S., 
for  which  Q  Is  lower  and  mor^  erratic,  and  the  eastern  U.S., 
for  which  Q  is  more  stable.  The  low  value  of  0  is  an  aver¬ 
age  value  and  may  result  from  both  low  and  high  Q  structures 
under  South  America. 
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A  Technique,  for  Determining  Earthquake  Source 
Parameters  by  Using  Surface  Wave  Amplitude  Spectra. 

by  Robert  B  Herrmann 

An  effective  method  for  the  determination  of  the  strike, 
dip,  and  slip  angles,  focal  depth,  and  seismic  moment  by 
using  the  fundamental  mode  amplitude  spectra  of  Love  and 
Rayleigh  waves  has  been  developed.  The  approach  taken  is 
to  find  the  particular  combination  of  these  source  param¬ 
eters  which  gives  the  best  agreement  between  the  observed 
and  theoretically  predicted  surface  wave  amplitude  spectra 

at  all  observing  stations  for  the  entire  period  range  of  the 
observed  spectra . 

Using  the  development  of  surface  wave  theory  by  Tsai 
and  Aki  (1970),  theoretical  Love  and  Rayleigh  amplitude 
spectra  can  be  generated  from  an  assumed  earth  model  for 
a  given  choice  of  the  five  source  parameters.  We  assume 
that  the  dislocation  source  is  a  step  function  in  time. 

^ eZ  yL^lj^  and  y^(TjP)  be  the  observed  Love  and  Ray¬ 
leigh  wave  amplitude  spectra  of  a  particular  mode  at  a 
period  T  and  azimuth  0  which  have  been  corrected  for  the 
geometrical  spreading  of  surface  waves  on  a  sphere  to  a 
reference  distance  of  o°  (i000  km)  and  also  corrected  for 


the  decrease  in  amplitude  due  to  the  effects  of  anelastic 
attenuation.  If  AQ  is  the  ooserved  amplitude  spectrum  of 
trie  Q  wave  type  at  a  distance  A  from  the  source  and 
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is  the  appropriate  anelastic  attenuation  coefficient, 
then  the  corrected  amplitude  spectrum  is 

yQ  =  [  sin  9"  V  exp{  !T«A)  • 

For  a  given  choice  of  the  source  parameters,  strike 
(STRK),  dip  (DIP),  slip  (SLIP),  and  depth  (H),  theoreti¬ 
cal  amplitude  spectra  can  be  generated  at  a  reference  dis¬ 
tance  of  9°  for  a  unit  dislocation  source  such  that  for 
each  observed  yq(T,j^),  a  corresponding  theoretical  amplitude 
spectrum  xq(T,0)  can  be  predicted.  Thus  for  each  data  point 
of  the  observed  amplitude  spectrum  specified  by  its  period, 
azimuth,  and  mode,  there  corresponds  a  predicted  amplitude 
for  the  same  period,  azimuth,  and  mode.  Having  established 
this  one-to-one  relationship,  we  will  drop  the(T,0)  descrip¬ 
tion  of  the  observed  and  predicted  spectra  and  just  consider 
the  yQj  and  Xqj  observational-theoretical  pair,  where  the  j 
Index  indicates  a  particular  yQ  -  Xq  pair  from  the  data  set 
of  all  observed  periods  and  azimuths. 

Let  there  be  a  total  of  NL  Love  wave  observations  of 
the  amplitude  spectrum  and  NR  Rayleigh  wave  observations. 

The  data  sets  consist  of  the  totality  of  Love  and  Rayleigh 
wave  spectral  amplitude  observations  available  in  the  range 
of  observed  periods  at  all  azimuths.  Using  the  simplified 
notation,  the  following  parameters  are  def.'ned: 
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MR  = 


NR 


NR 


E  *aj‘ 
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NR  .  MR  +  ML  »  MT. 
_NR  +  NL 


'NL 


RESL  =  JL  (yL1  -  no  ■  xT]y 

j=i  NL 


NR 
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E  <yH.1  -  m  • 
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-i  A 
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The  significance  of  these  parameters  is  as  follows. 
The  RL  and  HR  are  correlation  coefficients  between  the 
observed  and  theoretical  Love  and  Rayleigh  wave  data 

sets.  The  r.nrn«lnfln«  _  i j.  .  , 

- - **  — JL3  no  I,  affected  by 

multiplicative  scaling,  e.g.  seismic  moment,  but  rather 
gives  a  measure  of  the  similarity  of  the  observed  and 
theoretical  azimuthal  radiation  patterns  over  the  entire 
period  range  of  the  observations.  ML  and  MR  are,  re¬ 
spectively,  least  squares  estimates  of  the  seismic  moment 
using  the  Love  and  Rayleigh  wave  data  sets  independently. 
HO  is  the  weighted  estimate  of  the  seismic  moment. 


-  *-^1  "iirr  111  Ml  Hiiii 
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RESL  and  RESR  are  the  root  mean  square  amplitude  spectrum 
residuals  between  the  observed  amplitude  spectra  and  the 
theoretical  predictions  made  using  the  average  seismic 
moment  MO.  RES  is  the  root  mean  square  amplitude  spec¬ 
trum  residual  for  the  entire  data  set. 

The  technique  for  selecting  the  source  parameters 
which  best  describe  the  surface  wave  observations  follows. 

For  each  particular  choice  of  SLIP,  DIP,  STRK,  and  H,  the 
following  quantities  are  computed  and  displayed:  MO,  RR, 

RESR,  MR,  RL,  RESL,  ML  and  RES.  The  best  focal  mechanism 
and  focal  depth  combination  is  taken  to  be  the  one  for 
which  the  correlation  coefficients  RR  and  RL  are  as  large 
as  possible  and  for  which  the  seismic  moment  estimates  ML 
and  MR  are  also  as  close  as  possible.  Usually  it  happens 
that  for  the  source  parameters  satisfying  these  criteria, 

RESL,  RESR,  and  RES  are  also  minimal.  The  advantage  of 
this  algorithm  is  that  equal  weight  is  given  to  the  Ray¬ 
leigh  and  Love  wave  observations  and  that  noise  in  one  data 
set  does  not  unduly  affect  the  best  choice  made  by  the  other 
data  set.  A  disadvantage  is  that  the  systematic  search  is 
very  time  consuming  on  a  digital  computer,  as  opposed  to 
iterative  methods. 

The  search  is  usually  begun  over  the  following  param¬ 
eter  space: 


. ■< — - - - -  - - - - -  -  ■  -  . . -  -  - - - — 


SLIP 

-90° 

(20°) 

70( 

DIP 

u> 

o 

o 

(15°) 

90( 

STRK 

0° 

(20°) 

l6o' 

H 

0 

(  5  ) 

25 

The  numbers  within  the  parentheses  Indicate  the  Increment¬ 
ing  value  used  In  going  from  the  smallest  to  the  largest 
value  of  the  particular  parameter.  Due  to  the  symmetry  of 
the  surface  wave  radiation  pattern  to  a  l80°  rotation  for 
a  double-couple  seismic  source  and  because  only  amplitudes 
are  being  used,  these  search  limits  encompass  all  possible 
focal  mechanisms. 

Once  a  solution  region  is  located,  a  finer  search  grid 
is  used  until  all  angles  are  known  to  within  5°  and  the' 
depth  is  known  to  within  2  kilometers.  Usually  the  results 
of  this  technique  are  not  a  single  best  solution,  but  rather 
a  suite  of  possible  solutions  to  the  surface  wave  amplitude 
problem.  The  choice  of  the  focal  mechanism  of  the  earth¬ 
quake  is  taken  to  be  that  member  of  the  suite  of  possible 
solutions  which  gives  the  best  fit  to  the  P  wave  first 
motion  data. 


References 
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Appendix.  A  Summary  of  Central  United  States  Earthquake 
Source  Parameters  Determined  from  Surface 
Wave  and  P  Wave  Observations  . 

by  Robert  B.  Herrmann 

In  the  previous  section,  a  systematic  search  procedure 
for  the  determination  of  earthquake  source  parameters  using 
surface  wave  amplitude  spectra  was  discussed.  The  results 
of  applying  this  method  to  the  data  available  from  some 
central  United  States  earthquakes  are  presented. 

The  amount  of  data  available  for  the  surface  wave 
studies  depended  on  the  size  of  the  earthquake  studied  and 
the  number  of  seismograms  available  for  the  event.  The  14 
August  1965  event  was  recorded  at  eleven  long  period  seis¬ 
mograph  stations.  There  were  88  Rayleigh  wave  and  91  Love 
wave  spectral  amplitude  observations  available  in  the  2-12 
second  period  range.  On  the  other  hand  the  03  March  1963 
event  had  529  Rayleigh  wave  observations  in  the  7-50  second 
period  range  from  36  stations  and  393  Love  wave  observa¬ 
tions  in  the  8-50  second  period  range  from  28  stations. 

For  the  other  events  studied  there  were  on  the  average  250 
Love  and  Rayleigh  wave  spectral  amplitude  observations  in 
the  8-50  second  period  range  from  20  stations  . 

In  analyzing  the  data,  an  attenuation  model  which 
represented  an  average  of  the  models  of  Mitchell  (1973) 
and  Herrmann  (1973)  was  used  to  correct  for  the  effects 
of  anelastic  attenuation  on  the  observed  amplitudes.  The 


9 

seismic  source  was  assumed  to  be  of  the  step  dislocation 
in  time  type. 

The  results  of  our  investigation  are  presented  in 
Table  I.  This  table  lists  the  date  and  origin  time  of 
the  event  together  with  the  epicentral  coordinates.  The 
orientations  of  the  nodal  planes  and  the  orientations  of 
the  pressure  and  tension  axes  as  determined  ey  the  joint 
use  of  surface  wave  and  P  wave  data  are  also  listed.  The 
surface  wave  observations  made  possible  an  estimate  of  the 
focal  depth  and  the  seismic  moment.  The  mb  and  Ms  magni¬ 
tudes  for  each  event  are  from  Nuttli  and  Zollweg  (1974). 

As  a  result  of  the  investigations  made,  the  seismic  moment 
is  estimated  to  be  accurate  to  within  a  multiplicative 
factor  of  1.5,  the  focal  depth  to  within  2.5  kilometers, 
and  the  nodal  plane  strike  and  dip  to  within  10°. 

Further  studies  will  be  made  relating  the  seismic 
moment  and  magnitude  with  a  careful  appraisal  of  any  dis¬ 
crepancies  which  might  occur. 

References 
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2.  The  Seismic  Scaling  Law  and  Source  Mechanism 
Characteristics  for  Small  Magnitude  Events  In  the 
Central  United  States . 

by  Ronald  Street 


The  objective  of  this  report  is  to  outline  the  results 
of  recent  research  directed  towards  characterizing  earth¬ 
quake  activity  in  the  six-state  area  of  Illinois,  Missouri, 
Tennessee,  Arkansas,  and  Kentucky.  With  this  objective  in 
mind,  source  characteristics  of  some  eighty  earthquakes,  and 
the  focal  mechanisms  of  more  than  thirty  of  the  larger  of 
these,  have  been  determined  . 


The  prin  ipal  result  of  that  portion  of  the  study  con¬ 
cerned  with  source  characteristics  has  been  the  establish¬ 
ment  of  a  scaling  law  of  the  source-time  function  involving 
displacement  spectral  densities  from  the  excitation  of  the 
Lg  waves  (as  observed  on  short-period  records)  for  earth¬ 
quakes  whose  seismic  moment  ranges  from  101®  through  102/| 
dyne-cm.  The  principal  result  of  the  portion  involving  focal 
mechanisms  has  ’been  the  establishment  of  two  generalized 
trends  of  earthquake  mechanisms  which  appear  to  be  closely 
related  to  local  tectonic  features  such  as  the  Mississippi 
Embayment,  the  Ozark  Uplift,  and  the  Pascola  Arch.  Data 
for  both  studies  were  gathered  from  seismograms  of  the  Long 
Range  Seismic  Measurement  (LRSM)  stations,  the  World  Wide 


Standardized  Seismic  Network  (WWSSN),  Saint  Louis  University 
stations,  and  individual  stations  operated  by  nearby  univer¬ 


sities  (see  Fig.  l). 


Figure 
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Concerning  the  source- time  function,  the  procedure 
followed  was  first  to  determine  the  displacement  spectral 
density  of  an  earthquake's  coda  beginning  with  the  maxi¬ 
mum  trace  offset  after  the  P  phase  arrivals.  For  the  spec 
tral  analysis  codas  were  digitized  from  this  point  on  the 
trace  until  such  time  that  they  were  substantially  less 
than  the  initial  maximum  sustained  motion.  At  small  epi- 
central  distances  the  maximum  trace  offset  in  the  earth¬ 
quake's  coda,  excluding  the  P  phases,  was  the  Sg  phase, 
while  for  distances  of  approximately  100  km  or  more  the 
maximum  trace  offset,  and  the  point  at  which  digitizing 
was  begun,  was  the  Lg  phase.  However,  since  all  but  the 


very  smallest  of  events  (those  with  moment  <•  lo19  dyne-cm) 
were  principally  based  on  data  taken  at  cpicentral  distances 
greater  than  100  km,  the  derived  displacement  spectral  den¬ 
sity  will  commonly  be  referred  to  as  that  of  the  Lg  phase. 

Once  the  Lg  displacement  spectral  densities  were  de¬ 
termined,  the  observed  zero  frequency  levels  for  eight 
large  events  for  which  Herrmann  (1974)  had  determined  the 
seismic  moment  were  noted.  Then  using  as  the  definition 
of  seismic  moment 


J£: 

where  |  = 

l  £  = 


density 

shear  wave  velocity 
zero  frequency  level 
reference  distance  , 


a  new  parameter  designated  as  the  spectral  density  level 
(hereafter  denoted  by  S.D.L.)  was  empirically  equated  to 


a  suitable  modification  of  R.  Values 


the  seismic  moment  by 
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used  for  p  and  ^  were  2.5  gm/cm^  and  3.5  x  10^  cm  sec 
respectively.  Once  the  S.D.L.  had  been  corrected  to  agree 
with  the  seismic  moment  for  the  eight  larger  events,  the 
nature  of  R  as  a  function  of  the  eplcentral  distance  was 
Investigated  by  comparing  the  S.D.L.  from  stations  varying 
in  distances  of  20  to  400  km.  The  resulting  equations  of 
the  S.D.L.  for  events  studied  in  this  report  were  deter¬ 
mined  to  be 


S.D.L.  — 


1.34  (xlO25)  r/100 
1.34  (xlO25)  J  r/100 


r  •<_  100  km 


r  100  km 


(-Mo) 


where  r  is  the  eplcentral  distance  in  km  and  the  values 
of  ^  and  are  as  previously  indicated.  Physically  this 
means  that  the  displacement  spectral  density  determined  for 
an  event  is  corrected  as  a  body  wave  for  ep. central  dis¬ 
tances  less  than  100  km,  whereas  for  greater  distances,  it 
is  treated  in  the  manner  of  a  surface  wave  correction.  In 
this  manner,  the  gross  spectral  properties  of  more  than  eighty 
earthquakes  were  determined  so  as  to  be  consistent  with  those 
events  whose  seismic  moments  were  found  by  Herrmann.  Fig.  2 
illustrates  a  number  of  these  events.  Indicating  zero  fre¬ 
quency  levels,  corner  periods,  and  a  variety  of  falloff 
behaviors . 

p 

Fig.  3  is  a  plot  of  the  observed  w  spectral  corner 
period  vs  the  displacement  spectral  density  levels  for 
many  of  the  events  shown  in  Fig.  2.  In  the  case  of  an  event 
exhibiting  two  spectral  corners  the  w“^  falloff  was  extra¬ 
polated  to  the  zero  frequency  level.  The  slope  of  the  upper 
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and  lower  line  segments  in  Pig.  3  is  w“3 .  it  is  interest¬ 
ing  to  note  that  the  three  line  segments  indicated  are 
remarkably  similar  to  Aki's  (1972 j  revised  model  A  theoret¬ 
ical  scaling  law  (see  Aki's  Pig.  9).  However,  it  should 
also  be  noted  that  the  plot  obtained  from  observational 
data  for  this  report  differs  from  Aki's  in  an  important 
respect;  the  observational  plot  in  Figure  3  is  displaced 
towards  shorter  periods  and  at  a  greater  S.D.L.  when  com¬ 


pared  to  Aki's  model  A.  This  fact  is  significant  if  we 
are  to  relate  the  three  quantities  Ms,  m^,  and  M0;  two  of 
these  quantities,  M0  and  m^,  are  plotted  as  functions  of 
each  other  in  Pig.  4.  Also  shown  in  Pig.  4  is  the  variation 
of  mb  with  the  magnitude  of  the  same  event  for  the  0.33  and 
0.1  period  waves  in  addition  to  the  1.0  period  waves. 

The  next  step  was  to  determine  the  correctness  of  the 
magnitude  relationships  as  proposed  by  Nuttli  (1973)  by 
plotting  a  number  of  the  better  determined  events  as  a 
function  of  their  displacement  spectral  density  levels. 


Fig.  5  shews  the  result.  T 


ue  uxoxeu  anu  aOj.au  lines  are 


two  different  least  squares  fits  to  the  data  shown.  The 
broken  line  is  the  least  squares  fit  to  all  the  data,  while 
the  solid  line  is  the  least  squares  fit  to  the  data  consid¬ 
ered  as  being  made  up  of  two  groups  with  moments  greater 
and  less  than  10  dyne-cm,  respectively.  Using  either  of 

the  two  cases,  it  appears  that  for  events  whose  S.D.L.  is 

22 

less  than  approximately  10  dyne-cm,  the  magnitude  is  pro¬ 


gressively  overestimated  as  the  size  of  the  events  decreases. 


18 


■  a -  1  —  1  y 

DISPLACEMENT  SPECTRAL  DENSITY  (DYNE-CM) 


mb  vs  10  hz 


This  is  the  result  of  the  break  in  the  slope  of  the  m^  vs 
S.D.L.  scaling  near  the  1022  dyne-cm  S.D.L.,  and  the  fail¬ 
ure  of  Nuttli 's  magnitude  relationship  to  reflect  the  change 
in  the  displacement  spectral  density  spectrum  observec  at 

the  1.0  second  period  for  earthquakes  with  S.D.L.  less  than 
22 

10  dyne-cm.  This  overestimation  of  magnitude  could  sig¬ 
nificantly  modify  the  m^  vs  felt  area  (Nuttli  and  Zollweg, 
1973 )>  and  the  recurrence  relationships  (Nuttli,  1973)  for 
central  United  States . 

The  second  part  of  this  research  effort  involved  a  de¬ 
tailed  focal  mechanism  investigation  of  moderate  size  earth¬ 
quakes  occurring  during  the  last  thirteen  years  and  within 
the  six-state  region  of  interest.  The  events  studied  are 
listed  chronologically  in  Table  1.  Figure  6  illustrates 
the  focal  mechanisms  determined  for  these  events  from  the 
first  motion  of  P.  The  mechanism  solutions  shown  include 
those  of  Herrmann  (1974) ,  which  were  obtained  from  the  con¬ 
sideration  of  the  amplitude  spectrum  of  surface-wave  motion, 
as  well  as  the  thirty  solutions  based  on  the  sense  of  motion 
of  crustal  P  phases  (Pn,  P*,  Pg)  determined  for  this  report. 
For  those  events  based  on  the  sense  of  motion,  compressions 
or  dilations  were  plotted  on  the  lower  hemisphere  of  an 
equal  area  stereographic  projection  at  the  appropriate  azi¬ 
muth  and  incidence  angle.  Depending  on  the  particular  event 
it  is  felt  that  the  errors  in  the  strike  and  dip  of  the 
nodal  planes  are  no  more  than  10  -  20°. 
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Table  1.  Focal  mechanism  parameters  for  Central  United 

States  earthquakes. 


EVENT  DATE 


OT(UT) 


LAT 
(°  N] 


P  AXIS 

TREND  PLUNGE 


T  AXIS 
TREND  PLUNC 


02  FEB  62 

06  43  3^.0 

36.5 

89.6 

52 

26 

278 

01  JUN  62 

11  23  40.5 

35.0 

90.2 

336 

80 

252 

14  JUL  62 

02  23  49.0 

36.8 

89.9 

92 

2 

356 

23  JUL  62 

06  05  17.0 

36.1 

89.4 

131 

75 

233 

03  MAR  63 

17  30  11.4 

36.7 

90.1 

174 

10 

76 

31  MAR  63 

13  31  03.7 

36.9 

89.0 

161 

83 

104 

06  APR  63 

08  12  22.4 

36.5 

89.6 

85 

7 

346 

03  AUG  63 

02  37  47-8 

37.0 

88.7 

92 

3 

191 

16  JAN  64 

05  09  57.1 

36.8 

89.5 

220 

15 

95 

17  MAR  64 

02  15  29.0 

36.2 

89.6 

304 

32 

168 

23  MAY  64 

15  00  33-7 

36.6 

90.0 

223 

68 

83 

11  FEB  65 

03  40  24.0 

36.4 

89.7 

268 

10 

33 

06  MAR  65 

21  08  50.5 

37.5 

91.1 

260 

67 

165 

14  AUG  65 

13  13  56.6 

37.2 

89.3 

238 

27 

147 

21  OCT  65 

02  04  38.4 

37.5 

91.1 

340 

85 

160 

04  NOV  65 

07  43  33.6 

37.1 

91.0 

158 

5 

289 

12  FEB  66 

04  32  14.7 

35.9 

90.0 

140 

83 

271 

13  FEB  66 

23  19  36.9 

37.0 

91.0 

292 

68 

183 

26  FEB  66 

08  10  19.4 

37.1 

91.0 

292 

68. 

183 

06  DEC  66 

08  00  47.0 

38.9 

92.8 

277 

2 

10 

04  JUN  67 

16  14  13.6 

33.6 

90.9 

248 

7 

156 

21  JUL  67 

09  14  43.9 

37.5 

90.6 

314 

52 

50 

Table  1 

Table  1.  Cont'd, 


JT  DATE 

ot(ut) 

09  NOV  68 

17  01  42.0 

01  JAN  69 

23  35  36.2 

28  FEB  69 

13  10  13.1 

27  MAR  70 

03  44  29.5 

17  NOV  70 

°2  13  54.5 

12  FEB  71 

12  44  27.2 

13  APR  71 

T4  00  50.0 

18  OCT  71 

06  39  30.7 

01  FEB  72 

°5  42  10.0 

LAT  LON 
(°  «)  (°  W) 


38.0 

34.8 

37.9 
36.5 


88.5 

92.6 

88.6 

89.7 


29  MAR  72 
07  MAY  72 
09  JUN  72 
19  JUN  72 
12  JAN  73 
03  OCT  73 
09  OCT  73 


20  38  31.9 

°2  12  O8.5 

15  19.1 
05  46  14.7 
11  56  56.0 
°3  50  14.0 
20  18  26.8 


35.9 

37.7 
37.0 

37.9 

35.8 
36.5 


90.0 
9  0.4 

89.1 

90.5 

90.1 

89.6 


P  AXIS 

trend  plunge 


97 

330 

282 

264 


35.9 

90.2 

71 

38.5 

87.9 

78 

35.8 

90.1 

54 

36.7 

89.6 

108 

36.4 

90.8 

227 

36.2 

89.6 

180 

52 

269 

138 

70 

251 


1 

7 

2 

75 

67 

50 

72 

8 

60 

86 

72 

49 

22 

25 

13 


T  AXIS 

trend  plunge 


192 

228 

192 

15 

307 

344 

263 

12 

128 

270 

265 

15 

268 

278 

114 


82 

65 

78 

8 

13 

3 

15 

53 

10 

10 

15 

14 

60 

63 

72 
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Focal  mechanisms  done  by  Herrmann  (l,  5,  i4,  2.5,  21 , 
22,  23,  and  24)  were  obtained  by  the  assumption  of  a 
hypothetical  event  which  was  then  varied  until  the  theo¬ 
retical  radiation  patterns  for  Love  and  Rayleigh  waves 
compared  favorably  to  those  actually  observed.  A  con¬ 
straint  on  the  chosen  solution  is  that  it  must  be  in  rea¬ 
sonable  agreement  with  the  observed  sense  of  motion  of  the 
various  P  phases.  Errors  in  the  strike  and  dip  of  the 
nodal  planes,  as  determined  by  this  method,  are  probably 
less  than  10  .  While  this  is  obviously  a  superior  method 
to  that  which  uses  P  wave  data  only,  it  is  restricted  to 
events  which  excite  fundamental  mode  surface  waves  in  the 
4  to  50  second  period  range  sufficiently  so  that  they  are 
measurable  on  long-period  seismograms . 

The  strike  of  the  nodal  planes  shown  in  FJ.g.  6,  with 
a  few  exceptions,  can  be  grouped  into  two  trends.  The 
most  prominent  is  a  group  of  earthquakes  whose  nodal  planes 
strike  nearly  north-south  and  whose  epicenters  extend  north¬ 
ward  from  Memphis,  Tennessee  through  central  Illinois.  The 
second  grouping  of  earthquakes  whose  nodal  planes  appear  to 
correlate  with  one  another  are  those  in  central  Missouri. 

The  strike  of  the  nodal  planes  of  this  group  tend  northwest- 
southeast  . 

Two  significant  points  to  be  noted  about  these  two 
trends  are  the  complexity  of  their  intersection,  and  the 
change  In  character  of  the  north-south  trend  from  high- 
angle  normal  faulting  south  of  36.3°N  (near  the  buried 
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Pascola  Arch)  to  high-angle  thrust  faulting  to  the  north. 
The  focal  mechanisms  In  the  vicinity  of  the  Intersection 

of  the  two  trends,  as  indicated  by  the  mechanism  diagram 
on  the  figure,  is  complex. 

The  complexity  of  the  focal  mechanisms  near  the  Inter¬ 
section  of  the  two  trends  was  supported  by  an  attempt  to 
resolve  a  composite  fault  plane  solution  Involving  twenty- 
two  microearthquakes  occurring  in  the  Immediate  vicinity  of 
36.5°N/09./|°w.  biased  on  the  sense  of  first  motions  at  five 
seismographs  (ORV.  BRM,  DY3,  DY4,  and  DY5)  located  in  the 
area  surrounding  this  point,  it  proved  to  be  Impossible  to 
separate  the  composite  fault  plane  stereographic  projection 
into  any  reasonable  number  of  solutions.  This  supports  the 
contention  of  intersecting  fault  zones  as  well  as  implying 
something  other  than  a  simple  stress  field  as  the  cause  of 
observed  earthquake  activity  in  this  region. 


3  for  the  change  in  the  character  of  faulting 


north-south  trend,  the  implication  is  that  this  continuous 

Of*  r\v-\  4  4~  _____ 

" . ~^ccn„ers  is  chai'autez-iaed  b;y  a  ten- 

sional  stress  system  predominating  in  the  southern  portion 
and  a  compressive  stress  system  predominating  in  the  north. 
Three  mechanisms,  those  of  03  Mar  63,  01  Jan  69,  and 

14  AUS  65'  have  a  P03slbly  interesting  interpretation. 

All  three  occur  near  the  intersections  of  nodal-plane 
trends.  Pig.  7  illustrates  the  two  trends  previously 
described  and  the  three  strike-slip  focal  mechanisms 
Hated  above.  While  the  relationship  of  the  events  of 
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03  Mar  63  (No.  5)  and  14  Aug  65  (No.  14)  to  the  Intersec¬ 
tion  of  the  two  trends  is  clear,  that  of  the  01  Jan  69 
event  (No.  24)  is  not.  However,  if  Oetking's  (1966)  de¬ 
scription  of  the  Ouachita  Front  as  being  a  thrust  feature 
is  accepted,  then  the  mechanism  of  this  event  could  be  the 
result  of  a  superposition  of  the  thrust  fault,  compressive- 
stress  system  associated  with  the  Ouachita  Front  and  the 
perpendicular  normal  fault,  tensional-stress  as  found  at 
Memphis.  This  hints  at  the  possibility  of  tensional-type 
faulting  extending  southward  from  Memphis  along  the  axis 
of  the  Mississippi  Embayment  to  the  Ouachita  Front.  But 
at  present  the  distribution  of  seismograph  stations  is  not 
adequate  to  determine  if  earthquakes  are  occurring  along 
such  an  extension. 

The  results  of  the  research  presented  in  this  report 
bear  on  three  important  problems:  the  source  time  function 
of  micro-  and  moderate-size  earthquakes,  the  state  of  stress 
in  the  interior  of  a  continental  lithospheric  plate,  and  the 
relation  of  present-day  earthquake  activity  in  the  central 
United  States  to  geological  features.  With  regard  to  the 
first,  it  has  been  shown  that  it  is  possible  to  derive 
empirically  an  observational  scaling  law  which  is  consis¬ 
tent  and  which  is  directly  related  to  the  source  parameter 
of  seismic  moment.  This  capability  to  scale  properly  earth¬ 
quakes  in  central  United  States  will  have  a  significant 
impact  on  magnitude,  recurrence,  and  magnitude  vs  felt 
area  relationships,  and  indicates  the  need  of  establishing 
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similar  laws  appropriate  to  other  regions. 

With  regard  to  the  latter  two  problems.  It  has  been 

shown  that  the  stress  distribution  implied  by  the  observed 

focal  mechanism  solutions  is  considerably  more  complicated 

than  has  been  proposed  by  such  findings  as  those  of  Sbar 
and  Sykes  (1973). 
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3.  The  Q  Structure  above  the  Subduction  Zone 
in  South  America . 
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by  Rama  Somayajulu 


lateral  differences  in  the  anelastlcity  of  the  earth  is 
a  factor  which  affects  strongly  the  attenuation  of  seismic 
waves,  and  hence  the  magnitude  values  determined  from  seis¬ 
mic  amplitudes  recorded  at  seismograph  stations.  A  deter¬ 
mination  of  the  lateral  variations  in  the  Q  structure  of  the 
earth  is  thus  an  important  study  for  understanding  the  fac¬ 
tors  affecting  earthquake  magnitudes.  The  present  study 
is  an  investigation  of  the  anelastlcity  in  a  particular 
type  structure,  the  region  above  the  *:azca  Plate  in  the 

vicinity  of  the  Poru-Chlle  border  beneath  the' coast  of  South 
America . 

Tlx.re  are  two  classical  approaches  to  find  the  atten- 
uation  from  seismic  body  waves.  The  first  is  to  measure 
and  correlate  the  amplitudes  of  the  same  phase  from  station 

f-  r.  A.  J  „  . .  , 

lr‘e  aeU°nd  tu  ^ermine  the  amplitude  ratio 
Of  two  different  phases  from  the  same  source  recorded  at 

the  same  station.  The  former  method  involves  corrections 
roui  -e,  tne  .instrument,  the  receiver  crust,  and  the 
wave  paths.  In  the  latter  metnod  it  is  possible  to  ovoid 
tome  of  these  corrections,  namely,  those  for  the  instrument 
and  receiver  crust,  and  sometimes  for  portions  of  the  wave 
path.  Many  authors  have  applied  this  latter  method,  using 
the  spectral  ratios  of  the  direct  P  or  S  phase  and  the 
phase  transmitted  and/or  reflected  at  the  nor*  boundary 


fiSksiUii,, 
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or  at  the  surface  of  the  earth  (Kanamori,  1967;  Sacks, 

19o9;  Niazl,  1971;  Kovach  and  Anderson,  1964;  Steinhart 
et  al.,  1964;  etc . ) 

The  method  followed  In  this  study  is  to  examine  the 
spectral  ratio  of  successive  pp  and  P  phases  observed  at 
many  stations  from  the  same  source.  For  deep  focus  earth¬ 
quakes  and  for  epicentral  distances  greater  than  60°,  the 
direct  P  wave  travels  wholly  in  the  mantle,  whereas  the 
tirst  segments  of  pP  phase  traverse  the  crust  and  upper 
mantle  above  the  source  twice  before  travelling  a  path 
through  the  mantle  similar  to  the  direct  P  wave .  We 
assume  l)  that  the  spectral  difference  between  pP  and  P, 
after  correction  for  geometric  spreading  and  energy  trans- 
ler  at  the  surface,  is  produced  by  the  dissipation  properties 
of  the  crust  and  upper  mantle  above  the  source,  and  2)  that 
the  ray  P  and  that  portion  of  pP  below  the  depth  of  the 
focus  encounter  identical  Q  structure  in  the  mantle  and  at 
the  receiving  station. 


Method  of  Analysis 

The  method  of  analysis  is  described  by  Niazi  (1971). 
Briefly,  let  us  denote  the  spectral  amplitudes  of  P  and  pP 
by  Aj_  and  Ag,  respectively.  These  amplitudes  are  functions 
of  the  angle  of  incidence  i,  and  angular  frequency  w.  We 
assume  Q  Is  independent  of  frequency  (Anderson  and  Archam- 
beau,  1964;  Knopoff,  1964).  We  may  write  their  dependence 
on  the  various  source  and  path  elements  (Fig.  l)  as  follows. 


Ai  Ki,cf)  =  (1j  cf)  Gi(l)  G^l.w)  IiCi.w) 

a2  (w,l,cf)  =  S2  (l,q)  G2(i)  C2(i,w)  I2(i,w) 

R(l,w)  -kd 


where  S1?  s2  are  the  effect  of  radiation 

Gl>  G2  are  the  geometric  spreading 

G! ’  C2  are  the  effect  of  crust  and  upper  mantle 
excluding  the  portion  of  pp  (a CD,  Flg.l) 
path  above  the  focus. 

R  Is  the  reflection  coefficient  of  the 

shell  above  the  focus 

D  Is  the  distance  ACD  along  the  ray  above 

the  focus 

K  ls  average  attenuation  above  the 

focus 


K  = 


Q 

1 

T 

oL 


-l 


TT  /CX.  TQ 

specific  attenuation  coefficient 
is  the  angle  of  Incidence  at  the  source 
Is  the  bearing  of  the  ray  at  the  source 
is  the  travel  time  along  the  path  ACD 
the  longitudinal  wave  velocity. 

Since  for  deep  focus  earthquakes  recorded  at  larger  eplcen- 
tral  distances  p?  and  P  traverse  virtually  the  same  upper 
mantle  and  station  crust  (by  assumption  1),  and  are  recorded 


by  the  same  Instrument^ 


*1  =  I* 


C1  “  c2 
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Thus  amplitude  ratio  Is 


is  expressed 


fjg  =  s2  g2  r  Exp(-KD) 


*T 


^1G1 


-KD 


(3) 

/ 

The  g60metrlcal  «n  be  computed  for  a  particular 

ray  path  fro ■  the  following  expression  (sullen,  1963) 


G  = 


Sin  9 


r2Sin  A  cos  0 


de 


where  r  is  the  radius  of  the  earth 

A  eplcentral  distance 

6  angle  that  ray  leaves  the  source 

©o  the  angle  of  emergence  of  ray  at  the  earth's 
surface  . 

The  geometric  spreading  for  P  and  pP  are  small.  The  re¬ 
flection  coefficient  at  the  surface  is  computed  by  using 
the  Thomson- Haskell  matrix  method  (Haskell,  i96a).  A  com. 

Plate  computer  program  for  determining  the  pP  and  P  travel 
geometric  spreading,  and  energy  partition  has  been 
written  in  order  to  make  comparison  for  energy  which  has 
followed  pP  and  P  paths.  The  above  program  makes  use  of 
a  multilayered  system,  where  the  velocity  m  each  layer 

“  ?eSOrlbe<I  a  P°"'er  of  the  for*  V  .  arb  (Mohorovlclc 
!•«  with  a  and  b  constants  (Herrin  et  al.  1968).  Rewriting 
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J" 


equation  3  and  substituting  for  K  and  D 


lo„  .  ^  S1  G1 
°e  \  A1S2G2  R 


=  -KD 


where  Q  is  the  average  value  of  Q  for  the  crust  and  upper 
mantle  above  the  focus.  We  may  also  write  equation  (4)  ac 


A,  S,  G- ,  \ 

log*  -» &  - — -  \  =  bf 

*\£> 

where  _  r 


Hence,  tne  spectral  amplitude  ratio,  after  correcting  for 
the  geometric  spreading,  reflection  at  the  surface  and  the 
el  feet  of  radiation,  is  a  function  of  frequency  with  a 
negative  slope.  The  integral^  Jo  ^  ,  in  equations  (4) 

(u)  corresponds  to  the  time  taken  by  the  P  wave  to  reach 
the  reflecting  point  at  the  free  surface  and  to  return  to 
the  level  of  the  focus.  This  value  can  be  readily  read  from 
the  above  computer  program. 

More  precisely,  the  slope  of  the  graph  of  the  logarithm 
ol  the  corrected  spectral  ratio  versus  frequency  is  equal  to 
(Buchbinder,  1972) 


<3  I  * 

'  vr 
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where  the  Integral  path  pP+  is  the  path  A CD  (Pig.  l)  of  pP 
above  focus  pP-  in  the  path  DB  of  pP  belovr  focus,  and  P  is 
the  path  AB  of  P .  At  the  large  epicentral  distances 
( >  60  )  the  radii  of  penetrations  to  P  and  pF  differ  by 
approximately  250  kms,  where  both  the  rays  traverse  the 
lower  mantle.  By  the  assumption  2,  the  last  two  terms  in 

equation  (6)  will  be  approximately  equal  to  one  another. 
Then 

Q  *  Q 


Figure  1  Path  of  P  and  pP.  ACD  is  the  path 
of  Integration  for  pP+  . 
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Data 


Table  1 (see  also  Fig.  2)  lists  the  earthquakes  and 
recording  stations  used  in  this  study.  The  selection  of 
the  stations  was  based  on  the  azimuthal  distribution  and 
on  the  availability  and  clarity  of  the  records.  The  hypo 
are  located  in  the  zone  of  deep  focus  earthquakes 
in  Central  Chile.  This  zone  lies  beneath  the  familiar 
' Quiet-  zone  (250-550  kms  depth)  for  earthquakes  under 
Chile .  The  region  above  the  hypocenters  is  one  in  which 
the  dip  of  the  Benioff  zone  changes,  witn  a  possible 


splitting  of  the  descending  plate  into  two  or  more 
(Stauder,  1973). 


tongues 


The  data  used  are  taken  from  long-period  vertical 
seismograms  only  of  WWSSN  stations  distributed  across  the 
United  States .  The  onset  of  the  pP  phase  is  approximately 
two  minutes  after  the  beginning  of  the  P  phase  at  these 
stations  and  for  the  earthquakes  selected.  Amplitude  spec¬ 
tra  for  the  P  and  PP  phases  were  computed  using  a  past 
Fourier  Transform  program.  Since  the  time  difference  be¬ 
tween  P  and  pP  is  small,  we  can  assume  that  m  taking  the 
ratio  for  the  two  phases  the  noise  spectra  win  cancel  out. 

The  spectral  amplitudes  were  corrected  for  radiation 
effect,  geometric  spreading,  and  reflection  at  the  surface. 
The  corrected  spectral  ratios  (ApP/Ap)  were  plotted  against 
the  frequency  on  a  semi-logarithmic  paper.  Regression  lines 
were  then  determined  in  order  to  obtain  the  best  fit  throurh 


Table 


•s  •» 

• 

•» 

* 

cyfe 

p 

o  < 

< 

p  * 

<c 

p  m 

o 

1  ' 

<5  O 

Q 

Q  CO 

n 

•N  • 

•% 

•<  •» 

•* 

C5  <  CO 

S 

O  Oh 

s 

, 

DWW 

P 

P  O 

< 

Q  CO  3 

P 

Eh  CO 

< 

O 

•% 

•*  •% 

•% 

■H 

O  pq  0< 

C 

o  o 

p 

P 

P  P  O 

p 

CO  P 

o 

CO 

Eh  P  co 

< 

o  E 

u 

P 

CO 

•< 

•s 

cq  p  < 

U 

O  P 

*OiJ 

P 

CO  Eh 

pun 

Q 

PQ  «c 

C3  < 

<D 

-O 


00 

in 

00 

o\ 

in 

VO 

Lf> 

in 

to 

CO 


x: 

p 

e'¬ 

t'- 

oo 

a 

en 

C'- 

e'¬ 

a) 

Q 

in 

m 

en 

3 

o 


in 

in 

1 — 1 

cvj 

c 

• 

• 

• 

• 

o 

on 

on 

on 

on 

•H 

P 

CO 

vo 

VO 

VO 

vo 

« 

CO 

rH 

e- 

in 

O 

o 

• 

• 

• 

P 

CVJ 

s- 

t- 

e- 

OJ 

OJ 

CVJ 

CVJ 

<u 

ta 

E 

on 

vo 

rH 

•=t 

■H 

• 

• 

• 

• 

Eh 

rH 

f- 

-=f 

OJ 

in 

-=e 

■=t 

•3- 

C 

E 

•H 

vo 

e'- 

vo 

in 

bO 

on 

iH 

o 

on 

«H 

x: 

U 

OJ 

in 

o 

on 

O 

Ol 

iH 

rH 

rH 

CO 

OJ 

VO 

C- 

VO 

vo 

■3- 

•» 

vo 

•* 

Ov 

•» 

on 

OJ 

Ov 

•» 

CVJ 

0> 

CD 

p 

P 

P 

to 

a 

a 

o 

CO 

3 

<D 

0) 

a) 

Q 

< 

CO 

CO 

Q 

UW*---  -1 


38 


the  points.  The  slope,  b,  of  the  line  is  equal  to  the 

integral  term,  -  ^  |j>  ciS/^  J  inequation  (5).  Thus,  given 

the  value  of  b,  Q  can  be  calculated  for  each  earthquake  at 
each  station. 

The  calculated  mean  Q  values  so  obtained  are  listed 
in  Table  2.  Pig .  3  illustrates  these  same  values  together 
with  error  bars  indicating  the  standard  deviations. 

Event  I  has  Q  values  which  vary  from  100  to  580. 

Events  II  and  III  yield  Q  values  with  slightly  larger  varia¬ 
tions,  from  16  to  600.  For  event  IV  the  values  are  fairly 
uniform,  ranging  only  from  200  to  360.  Finally,  for  event  V 
Q  takes  on  values  notably  higher  than  in  the  other  cases, 
ranging  from  560  to  1440. 

In  Table  2  and  in  Figure  .3  the  receiving  statj  ons  are 
ranged  from  west  to  east,  numbers  1-6  corresponding  to  those 
in  the  western  United  States,  numbers  7-17  to  those  in  the 
eastern  United  States.  Examination  either  of  the  table  or 
the  figure  shows  that  the  Q  values  determined  from  stations 
in  the  western  United  States  are  both  lower,  on  the  average, 
and  more  erratic  than  those  determined  from  stations  in 
the  eastern  United  States.  While  the  reflection  point  of 
pP  and  the  precise  ray  path  of  Pp  relative  to  the  position 
of  the  slab  are  not  here  presented,  reference  to  Fig.  1 
Indicates  that  the  p  segment  of  pP  paths  for  the  western 
stations  passes  through  a  zone  just  above  the  upper  surface 
of  the  descending  oceanic  plate,  whereas  for  stations  to 
the  east  this  segment  of  the  pP  path  travels  through  greater 


17  SB fi  411.2  631.0  304.8 


portions  of  the  wedge  (aseismic  in  this  portion  of  South 
America)  between  the  descending  plate  (or  segment  of  palec 
plate)  and  the  continental  plate.  Evidence  from  seismicity 
and  focal  mechanisms  (Stauder,  1973)  shows  that  the  former 
path  is  likely  to  be  one  of  greater  heterogeneity,  for  the 
plate  is  probably  segmented  and  irregular  in  its  relief, 
thus  accounting  for  the  lower  ^  values  and  the  scatter  in 
the  values . 

The  values  of  Q  from  stations  in  eastern  United  States, 
while  more  uniform,  are  relatively  low  =  300  to  400) . 

This  indicates  passage  through  a  low  velocity  layer  of  con¬ 
siderable  vertical  extent  from  the  source  region  of  600  km 
depth  to  the  lithosphere.  Sacks  and  Okada  (1973 )  found  Qp 
values  in  excess  of  1000  generally  to  prevail  in  the  wedge, 
above  the  descending  Nazca  plate  at  depths  down  to  350  km, 
but  low  values  of  Qp  below  this  level .  The  values  reported 
here  are  average  values  for  the  entire  region  above  the  source 
level  at  600  km  depth,  and  were  obtained  from  a  double  passage 
(pP  +  path  above)  through  the  portion  of  mantle  in  question. 
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